Electronic Bloch oscillation in a pristine monolayer graphene 
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In a pristine monolayer graphene subjected to a constant electric field along the layer, the Bloch 
oscillation of an electron is studied in a simple and efficient way. By using the electronic dispersion 
relation, the dynamic formula of semi-classical velocity is derived analytically, and then many aspects 
of conventional Bloch oscillation, such as the frequency, amplitude, as well as the direction of the 
oscillation, are investigated. It is interesting to find that the electric field should affect the component 
of motion which is non-collinear with electric field, making the particle accelerated or oscillated in 
another component. 



Since its discovery in 2004, graphene has attracted a 
tremendous amount of interest due to its unique prop- 
erties that may promise a broad range of potential 
applications |l|-|5(. Recently, many theoretical and ex- 
perimental investigations focus on the graphene-based 
supcrlattices with electrostatic potentials or magnetic 
barriers 0" 11 1, including periodic [l2 - 1 5 1 . aperiodic [TgI— 
n~8j | . and disorder system Some researchers have also 
investigated the electronic Bloch oscillations in a struc- 
ture with periodic potentials fioj and a graphene nanorib- 
bon in presence of a periodic mass potential [2l| . One 
important issue still remains opening is that, what does 
the electronic Bloch oscillation behave in the pristine 
graphene without any potential structures or nanorib- 
bons? 

Bloch oscillation is an important phenomenon. It is 
usually involved with the coherent motion of quantum 
particles in periodic structures. For example, an electron 
(a matter wave) suffers this effect in a periodic lattice 
subjected to a constant external field. This phenomenon 
is predicted from quantum mechanics in very early 
days[22|, and has been demonstrated in various fields 



of physics, such as semiconductor super lattices [2J, |25j . 
photonic crystals 27 1, cold-atom systems [2^, and acous- 
tic waves [29j. Different from the common semiconduc- 
tors, the electron in graphene are described by Dirac 
rather than the Schrodinger equation, and then its Bloch 
oscillation naturally is an interesting issue to be investi- 
gated. 

Many aspects of Bloch oscillation can be obtained by 
a single band description via using the dispersion rela- 
tion to derive the semi-classical velocity of the particle. 
In this paper, based on the electronic structure under 
tight-binding approximation, we derive the motion of an 
electron in pristine monolayer graphene subjected to a 
constant external field. Within such a simple and effi- 
cient way, our results show several interesting phenom- 
ena of the electronic Bloch oscillation in graphene. In the 



following, we firstly derive the general formula of the mo- 
tion of an electron based on the dispersion relations, and 
then we analyze the properties of the Bloch oscillation. 

A monolayer graphene is well known for its honeycomb 
structure, and its dispersion relation can be written asQ 



8{k) = ±ey/3 + f(k), 



(1) 



2cos(v / 3afc y ) + 4cos(^-ak y )cos(^ak x ). 



where f\™) — i luo^ y ULin^ -r — ^~ u>nj y j 

a w 1.42 A is the carbon-carbon distance, and e ~ 
2.5eV is related to Fermi velocity (vf ~ 10 6 m/s), 
Hvf = feaQ- The signs "+" and "— " are, respec- 
tively, corresponding to the electron and hole energy 
band, which touch together at Dirac Points (DPs) . From 
Eq. [T] it is easy to find that the DPs are located at 
r],[^(2n+l)7r,^(2n±j)7r] with 

we can 



f^,-^(2n± i) 

1 3a ' ^/3a y 3 > 
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readily have v = (v x , v y ) as the function of k x and k 



According to v = ^ d£ S® t3C 



dk 



<j- 



=F3ea sin(|afc 2; ) cos(^a/cy) 
= FV / 3£a[sin(v / 3afc i) ) + sin(^afcj,) cos(|afc x )] 



hy/3+f(k) 
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which show that v x and v y 
of k x and k y , and the sign 



are the periodic functions 
'— " ("+") is corresponding 



to the velocities of electron ( hole or hole-like electron). 
When a constant electric field E = (E x ,E y ) is applied 
along the layer of graphene, the electronic motion equa- 



tion h dk jP = —eE survives as 



dt 



eE 

k x (t) = MO) ~ -j^t, 



k y {t) 



k y (0)-^t, (3) 



where £^(0) and k y (0) are the initial wave- vector values. 
Substituting Eq. [3]into Eq. 2, we can obtain the dynamic 
formula for v. Therefore using Eqs. 2-3, we can analyze 
the motion of the electron or hole in graphene under the 
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constant electric field. Since the direction of the electric 
field E can be chosen arbitrarily, we shall firstly discuss 
the electronic motion when E is only along the x (case 
I ) or y (case II) direction and then generalize it to an 
arbitrary direction (case III). 

Case I: E along the x direction. In this case E y = 0, so 
we have k y (t) = k y = constant, and k x (t) = k x (0) — ^j^t. 



Assuming k x (0) 
are 



0, the dynamic formula of v x and v y 



v x (t) 
v y (t) 



=F3ea sin(-^t) cos(^-ak y ) 



hG(t) 

=fy/3£a[sm(y/3ak y ) + sm(^aky) cos(— %ft)] 
h~G(t) ™ 



where G(t)=J 3 + 2 cos( % /3afc. y ) + 4cos(-^i) cos(^a/c y ), 
and T = i ^ , . From Eq. 4, it is easy to see that 
v(t + T) = v(t) with T being the period of the motion. 
The frequency and circular frequency of the Bloch 
oscillation are generally given by 



1 
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(5) 



respectively. According to the expression of v, the time- 
dependent position r(t) of the electron is r(t) = r(0) + 
J vdt, and here we assume the initial position r(0) = 0, 
i.e., x(0) = and y(0) = 0. After a simple derivation, 



we obtain x(t) = C 



tE 



■G(t) where C is an integra- 



tion constant satisfying xfo) = 0. For y(t), we have to 
numerically calculate the following 



y(t) = 



2^/Seujb 
3eE x 

* sin(v / 3afe y ) + sm(^a£;y) cos(— ujst) 
/o W) 
According to the formula of x(t), we can have 



eft (6) 
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eE : 
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Therefore, the amplitude of the oscillation along x direc- 
tion, L x = |x max - x mi „|, is given by 



,V3 



,V3 



Lx = |^ I|1 + 2C ° S( ~ afc »)|-| 1 - 2cos (-2- ak v)W- ( g ) 
When | cos(^fc. y a)| > i L x has its maximum value, 
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7§-t\ and when ^ak y = (n + \)tt, L 



0. 



According to Eq. [SJ if E x = 4.61 mV/nm[2l|, ( we set 
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FIG. 1: (Color online) The time dependence of Vx (a), t> y (b) 
and (c) the trajectories of the electron (or hole-like electron) 
on the graphene sheet for different k y at ^ak y = ( solid 
dark line line ), ^ak y — j (dash-dot-dot magenta line), 
^ak y — ^ (dash red line), ^ak y = -^ir (dot green line), 

and ^ak y = ^ (solid blue line). The electric field E is along 
the x direction. 



e = 2.5 cV in the whole paper ), £™ ax ~ 1084 nm with 
v B ~ 237GHz. 

Fig. 1(a) and (b) demonstrate the time dependence of 
v x and v y with different values of k y . It is found that 



when ^aky 



■|, shown as the dash red lines, the elec- 



tron shall pass through the DPs, and the period of v x and 
v y shall be doubled. Because the electron passes through 
the DPs, the electron should transit into another band 
and behave as a hole-like electron. After a period in an- 
other band, the hole- like electron shall behave as the elec- 
tron again. Therefore the period of the velocity becomes 
twice time, and correspondingly the amplitude shall be 
doubled. There is an interesting phenomenon that, the 
oscillation along the x direction disappears although E 
is still along the x direction when ^ak y = (n + \)k. 
Meanwhile, since v y ^ 0, the oscillation in the y direc- 
tion remains, see the solid blue lines (^ak y = ^) in 

Fig. 1(a) and (b). In other hands, if ^ak y = nn, we 
shall have v y = 0. It means that the oscillation in y-axis 
disappears and the oscillation in x-axis remains, see the 
solid dark lines in Fig. 1(a) and (b). 

The corresponding electron's trajectories arc shown in 
Fig. 1(c), where we demonstrate the trajectory of an 
electron within three periods on the graphene layer. It 
is clear that, when the electron (or hole) passes through 
the DPs, its amplitude shall be doubled and its trajec- 
tory is approximately a circle, see the dash red lines 
(^■ak y = -|) in Fig. 1(c). In general, the motion in 
the x and y directions may oscillate, but its trajectory 



3 



is very complex and depends on the initial value of k y . 
For example, its trajectory should be a helix (see short 
dash green line for ^ak y = y^7r); or go further and 
further with variational velocity, its trajectory should be 
like a sine function (see dash-dot-dot magenta line for 

2 ul% v ~ 4J- 

Case II: E along the y direction. In this case, E x = 0, 
so we have k x (t) = k x = constant and k y (t) = k y (0) — 
^"■t. Assuming k y (0) = 0, the dynamic formula read as 



V x (t) 

Vy(t) 



=F3ea sin ^(ak x ) cos(-^-t) 

HG'{t) 
Tv / 3ea[sin(^t) 



sin(-^Lt) cos(^ak x )] 



HG'(t) 



(9) 



where G'(t) = ^3 + 2 cos(-fR) 4 
and the period of the motion T' - 
M we still have v(t + T') = v(t) 



Acos(^ak x ) cos(-|ft) 
From Eq. 
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The frequency and 
circular frequency of Bloch oscillation are accordingly 
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and uj'-r = 
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FIG. 2: (Color online) The Lissajous figures of v x (a) and 
v y (b) for different values of k x , and the trajectories of the 
electron (or hole-like eletron) on the graphene sheet are shown 
in (c) and (d) with ^ak x = (solid dark line), ^ak x = | 
(dash magenta line), ^ak x — -| (solid red line), ^ak x = 

j^-k (dot green line), and ^ak x = ^ (solid blue line). The 
electric field E is along the y direction. 



From Eqs. 5 and 9, we can see that the direction of the 
electric field shall affect the frequency of Bloch oscilla- 
tion. Similar to Case I, y(t) = C - jf-G'(t), where C 
is an integration constant satisfying y(0) = 0. For x(t) 
we have to numerically calculate 



x(t) 



2^/3euj' b f* sin (^ak x ) cos {—u)' B t) 



eE v 



G'(t) 



dt (11) 



The amplitude of the oscillation along the y direction, 
L y = |y max - Vmia\, is given by 



Li, 



\eK, 



\4:Cos(-ak x 



cos 2 (-ak x ) 



(12) 



which show that L y has its maximum value L™ ax 



3e 



when cos(|afc x ) 



±1, and L™ m 



0. If E v 



mm 



671 nm with v' B 



4.61mV/nm, then L™ ax w 
137 GHz. 



\eE y \ 

cos(|afca;) 
1627 nm, L\ 

Fig. 2(a) and (b) show the Lissajous figures of v x and 
v y with different values of k x . When ^ak x = 0, the 
electron shall pass through the DPs and shall behave as 
a hole-like electron, and after a period of time it shall be- 
have as an electron again, and its amplitude is increased 
accordingly. At this case when — u't G [2n+|7r, 2n+|7r], 
the particle behaves as a hole-like electron, otherwise it 
behaves as an electron. 

Different from the case I, the oscillation along the y 
direction never disappears when electric field is along the 
y direction, and its amplitude never becomes zero. There 
is a special case that when k x = 0, we find that v x = 



0. At this case the oscillation in the ^-direction shall 
disappear. The corresponding trajectories within three 
periods are shown in Fig. 2(c) and (d). 

Case III: E along the arbitrary direction. At this case, 
the dynamics of the electron becomes much more com- 
plicated, since both E x and E y are not 0. Meanwhile, 
the dynamic properties of v x (t) and v y (t) are also re- 
lated with the initial phase k x (0) and k y (0) and the ratio 

Ex/Ey 

v x and v y depend on the two periodic functions 



\a[k x (0) 



< E 

n 



■t}) and cos(^a[fc !/ (0) 



tE 



^t] 



Let 



T x denote the period of the former, and T y denote the 



period of the latter, so we have T x = ^ i 



and 



rp 4lT H 

y ~ V3\aeE y \- 



If the ratio T x /T y is rational, i.e. 
Tx/Ty = m/n, (n and m are integers), then v x {t) and 
v y (t) are periodic with the periods being mT y or nT x . 
But if T x /T y is irrational, v x (t) and v y {t) should not be 
periodic anymore: they should not have a finite period. 
In this case the motion is not a periodic oscillation, even 
though the particle may move back and forth. 

Finally, we discuss the dynamics of electrons under the 
condition of k x (0) = k y (0) = with some specific di- 
rections, where a;(0) = y(0) = 0, a = arctan(i? y /i? x ), 
E = \E\, and lo x = 2tt/T x . The time dependence of v x 
and v y with different values of a are illustrated in Fig. 
3(a) and (b), and the corresponding trajectories within 
3Ta; are shown in Fig. 3 (c). Due to the rotational sym- 
metrical structure of graphene, it is easy to find that, 
when a = -| (solid dark line) the direction of the elec- 
tric field is equivalent to the y direction, and the electron 
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FIG. 3: (Color online) The time dependence of (a) v x , (b) v y 
and (c) the trajectories of the electron (or hole-like electron) 
on the graphene sheet for a = -| (solid dark line), a = \ (dash 
blue line) and a = ? (solid red line) . 



shall also pass through DPs in this case; when a = | 
(solid red line) the direction of the electric field is equiv- 
alent to the x direction. When ^ak y — j (dash blue 
line), it is a general case, in which the electron is only 
moving within a single energy band. 

In summary, we have derived the general formulas for 
the velocity of the electron in graphene subject to a con- 
stant electric field, and analyzed the dynamic proper- 
ties of electron for some particular and interesting cases. 
When electric field is along x-axis and y-axis, we find 
Bloch oscillation in direction of electric field, and we ad- 
dress formulas for its amplitude and frequency We also 
find that the electric field affect the motion in other direc- 
tion, making the electron oscillating or moving forward 
with fluctuation in other direction. Moreover, the veloc- 
ity is periodic in all directions in these two cases. Finally, 
we analyze the period of the motion and give the numer- 
ical result if electric field has an arbitrary direction. Our 
result should be a positive insight for experimentally ob- 
serving the Bloch oscillation in a pristine graphene, which 
may facilitate the development of graphene-based elec- 
tronics. 
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